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Abstract 

We present a model for neutrino oscillations in the presence of a deconstructed non- 
gravitational large extra dimension compactified on the boundary of a two-dimensional disk. 
In the deconstructed phase, sub-mm lattice spacings are generated from the hierarchy of 
energy scales between ~ 1 TeV and the usual B — L breaking scale ~ 10 15 GeV. Here, small 
short- distance cutoffs down to ~ 1 eV can be motivated by the strong coupling behavior 
of gravity in local discrete extra dimensions. This could make it possible to probe the 
discretization of extra dimensions and non-trivial field configurations in theory spaces which 
have only a few sites, i.e., for coarse latticizations. Thus, the model has relevance to present 
and future precision neutrino oscillation experiments. 
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1 Introduction 



There may be many ways to give a sensible short- distance definition of a model, which de- 
scribes the Kaluza-Klein (KK) modes [1] of a compactified higher-dimensional gauge theory. 
One attractive possibility is offered by deconstructed [2] or latticized [3] extra dimensions. 
Deconstruction provides a new type of four-dimensional (4D) manifestly gauge-invariant 
and renormalizable field theories^, which simulates in the infrared (IR) the physics of extra 
dimensions, and thus, yields a possible ultraviolet (UV) completion of higher-dimensional 
gauge theories [5]. In deconstruction, a bulk gauge symmetry is mapped onto iV copies of a 
4D gauge group with bi-fundamental Higgs "link" fields connecting the neighboring groups. 
This becomes equivalent with treating the extra dimensional space as a transverse lattice [6], 
where the inverse lattice spacing A is set by the vacuum expectation values (VEV's) of the 
link fields. These models find a graphical interpretation in "theory space" [7], where the 
notion of extra dimensions has been translated into more general organizing principles of 
4D field theories [8]. With this emphasis, the idea of theory space has stimulated the devel- 
opment of interesting theories in four dimensions, which need not be attributed any direct 
extra-dimensional correspondence at all. 

Since deconstruction is formulated within conventional 4D field theory, one may now 
expect that a typical lattice cutoff A naturally varies in the UV desert between ~ 1 TeV and 
Mpi ~ 10 18 GeV, i.e., in the energy range where usual effective field theories are defined. 
Indeed, if A is, e.g., of the order the B — L breaking scale Mb-l — 10 15 GeV, then a realistic 
neutrino phenomenology can arise [9], where small neutrino masses and bilarge leptonic 
mixing result from the lattice-link structure of the theory itselfH Generally, deconstructed 
non-gravitational extra dimensions always give a sensible effective field theory up to energy 
scales of the order 47rA, and hence, it seems that low cutoff scales A <C 1 TeV are not 
particularly preferred in this case. However, for inverse lattice spacings above a TeV it will 
be difficult to test theory space models at low energies, unless the number of lattice sites N 
is very large. On the other hand, cutoff scales which are hierarchically small compared to 
~ 1 TeV emerge in local theory space formulations of gravity. Local discrete gravitational 
extra dimensions are characterized by an intrinsic maximal inverse lattice spacing A max , 
which is related by a "UV/IR connection" to the higher- dimensional Planck scale M* and 
the compactification scale R^ 1 [11-13]. Owing to this UV/IR connection, one obtains a 
small cutoff A max <C M*, when R' 1 and M* are clearly separated. Low string scale models 
with M* ~ 1 TeV, like the large extra dimensional scenario of Arkani-Hamed, Dimopoulos, 
and Dvali (ADD) [14] (see, alternatively, Refs. [15,16]), will thus lead to a hierarchy A max <C 
1 TeV, when formulated in local theory spaces. 

This could, in principle, open up the possibility to test these models experimentally al- 
ready at low energies. For instance, it is well known that sub-mm sized extra dimensions 
would allow for a measurable conversion of the active neutrinos in higher- dimensional neu- 
trino oscillations [17-22]. Therefore, if discrete gravitational or non-gravitational large extra 
dimensions exhibit a cutoff A max ~ i? -1 ~ (mm) -1 , one could study with neutrino oscilla- 
tions the finite discretization effects of few site models far away from the continuum limit. 

1 For an early formulation in terms of infinite arrays of gauge theories, see Ref. [4]. 

2 A recent attempt to obtain the quark mixing from deconstruction has been described in Ref. [10]. 
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Clearly, in any such model, it would be important to have a dynamical understanding of the 
smallness of the inverse lattice spacing A < A max in terms of a mechanism, which generates 
a value A ~ (/im) -1 from energy scales in the UV desert above a TeV. 

In this paper, we use neutrino oscillations to probe deconstructed non-gravitational large 
extra dimensions, which have inverse lattice spacings of the order (/xm) _1 . We motivate 
the relevance of small inverse lattice spacings <C ITeV by the strong-coupling behavior of 
gravity in a six-dimensional (6D) model, where the two extra dimensions have been naively 
discretized in a local theory space. In fact, by requiring that a model of this type with only 
a few sites allows for testable predictions in typical Cavendish-like (laboratory) experiments, 
the maximal strong coupling scale A max is found to be A max < 1 keV. We analyze a toy 
model with few sites for a deconstructed U(l) gauge theory on a disk, where a sub- mm 
sized boundary emerges dynamically from the hierarchy of energy scales between ~ 1 TeV 
and Mb-l — 10 15 GeV. The active Standard Model (SM) neutrinos mix with a latticized 
right-handed {i.e., SM singlet) neutrino, which propagates on the boundary of the disk. 
A gaugeable cyclic discrete symmetry ensures that the latticized neutrino can be treated 
as a massless Wilson fermion. This symmetry also establishes a Wilson-line type effective 
coupling between the active neutrinos and the latticized right-handed neutrino. Hence, 
the model reproduces for coarse latticizations major features of the 5D ADD continuum 
theory with one right-handed neutrino in the bulk, which couples to the SM through a 
local interaction. We study the possible neutrino oscillation patterns for the cases of a 
twisted/untwisted right-handed neutrino and an even/odd number of lattice sites. The 
neutrino oscillations may account for possible subdominant deviations of neutrino oscillations 
from the solution to the solar and atmospheric neutrino anomalies and could be tested in 
ongoing and future low-energy neutrino experiments. 

The paper is organized as follows. In Sec. [H we review the strong coupling behavior of 
gravity in local theory spaces for the example two discretized gravitational extra dimensions. 
Next, in Sec. [3l we present our model for the deconstruction of a U(l) gauge theory on the 
boundary of a two-dimensional disk. This model generates from the hierarchy of energy scales 
between ~ 1 TeV and Mb-l lattice spacings in the sub-mm range, and thus, represents a few 
site model for large extra dimensions. Properties of the mass spectrum and the coupling of a 
latticized right-handed neutrino which propagates on the boundary of the disk are analyzed 
in Sec. HJ The mixing of the KK modes of the latticized neutrino with the active neutrinos 
is calculated in Sec. [5j In Sec. [61 we consider the neutrino oscillations of the active neutrinos 
into the latticized fermion and discuss the different neutrino oscillation patterns. Finally, in 
Sec. 0, we present our summary and conclusions. In Appendix |A] we study the cancellation 
of anomalies in our model, and in Appendix [BJ we describe in some more detail the neutrino 
mixing matrices. 

2 Strong coupling in discretized 6D gravity 

In this section, we will estimate the strong coupling scale in a naive discretization of 6D 
gravity. It turns out that in local theory spaces with sub-mm sized extra dimensions, the 
short-distance cutoff, as determined from graviton scattering, lies in the sub-MeV-range, 
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which is far away from the (effective) 6D Planck scale Mqd ~ 1 TeV. If we require that 
such a model leads to predictions which are testable in Cavendish-like experiments, then the 
UV cutoff can be even further lowered by several orders of magnitude. This motivates to 
consider the extreme limit of lattice cutoffs in the eV-range. This section briefly reviews part 
of the formulation of gravity in theory space as given in Ref. [11] by considering the case of 
two discrete gravitational extra dimensions. In doing so, we follow closely the treatment of 
a single discrete gravitational extra dimension in Refs. [12,13]. 

The starting point of our discussion is standard general relativity in six dimensions, 
where the two extra spatial dimensions have been compactified on a square. We write the 
coordinates in the 6D space as zm = (x^yk), where the 6D Lorentz indices are denoted by 
upper case Roman letters M = 0,1, 2, 3, 5, 6, while we use for the usual 4D Lorentz indices 
lower case Greek letters fi = 0,1,2,3, and the coordinates (k = 1,2) describe the fifth 
and sixth dimension. We write the 6D metric GMN{%fj,,yk) in block form as 

where we have neglected the associated spin-1 and spin-0 excitations^]. For our purposes, it is 
sufficient to restrict Gmn to the simplified form shown in Eq. (pQ), since we will be only con- 
cerned here with the leading order UV behavior of the scattering amplitudes in the effective 
theory for massive gravitons. The 6D Einstein-Hilbert action 5 EH = / d e x M^ D ^J\G\R & £,[G\, 
in which R n £,[G] is the Ricci scalar in n dimensions, can then be rewritten as 

Seh = f d A x d Vl dy 2 y^\MU^D[g] + -dy.g^g" 13 - g^g^)d yk g af s), (2) 

where summation over k — 1, 2 is understood. In order to simulate the effects of the two 
extra dimensions in a 4D model, we assume N 2 copies of 4D general coordinate invariance 
(GC), which we denote as GCi, where i = i 2 ) and i±, i 2 — 1, 2, . . . , N. In theory space, 
each coordinate invariance GCi is represented by one "site" i, where any two neighboring 
sites GCi and GCj with \j — i\ — 1 are connected by a link field Y^. Thus, the collection of 
sites and links forms a two-dimensional transverse lattice (see Fig. [1]). As already pointed 
out in Ref. [3] , this can be simply interpreted as a variation of the Eguchi-Kawai plaquette 
model for large N gauge theories [24]. 

Each site i is equipped with its own metric g^ixi), which is a function of the 4D coor- 
dinates Xi = (xf ) for this site. In order to obtain a lattice version of the derivatives d Vk in 
Eq. (T5]), we will employ appropriate pullbacks of the coordinates and metrics g l as defined 
in Ref. [11]. In this line of thought, a link field Yji can be considered as a pullback function 
Y^Xi), which maps a vector Xi on the site i onto a vector Xj on the site j with coordinates 
Xj = Y^Xi). Likewise, to compare the metric g^i^Xj) on the site j with the metric g^i^Xi) 
on the site i, one can introduce the field G%{xi) = {dY^/dx^dYfjdx^^^Xi)), which 
transforms as a metric tensor under GC,, but is left fixed by GCj. Therefore, G^ 1 is the 
pullback of the metric gi^Xj) on the site j to the site i. 

3 The spin-1 states decouple from the fields confined to the 3-brane, while the spin-0 fields interact only 
through the dilaton mode [23] . 
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Figure 1 : Plaquette in the theory space for a naive discretization of 6D gravity. Each circle denotes 
one general coordinate invariance (GC). Two neighboring GC's GC a and GCb, where \b — a\ = 1, 
are connected by a link field Yj, a , which can be regarded as a map from site a to site b. 

Assuming equal lattice spacings a in both the y\- and ^-directions, we can now associate 
the derivatives in the extra dimensional continuum theory with the nearest neighbor forward 
difference operators dg^j dy k = or 1 [g^fa) - G j * v (xi)] , where (ji, j 2 ) = (ii + 1, 12) for k = 1 
and (ji, j'2) = (*i,^2 + 1) for k = 2. With this definition, we obtain from Eq. ([2]) a naive 
transverse lattice formulation of the 6D Einstein-Hilbert action 

S = I d%^\M 2 {RM9 L ]+a~\g l ,M)-GU^)) 

x (g^( Xl ) - G^x.W^g^ - r a g^)) , (3) 

where \j — i\ = 1 and M is the fundamental 4D Planck scale, which can, in general, be 
different from the usual 4D Planck scale Mp\ ~ 10 18 GeV. Comparison with Eq. ([2]) shows 
that R = Na, where ~ 1/R is the compactification scale, Mg D = M 2 a~ 2 , and Mpi = NM. 
From Eq. ([3]) it is evident that we can now equally well drop on all the coordinates Xi the 
index i by making the identification = x^ and describe on all sites i the positions using 
only a single universal coordinate system with coordinates x M . 

Following the effective field theory approach advanced in Refs. [11-13], we can parame- 
terize the link fields Y£ in terms of small deviations from x M as 

Yjt(x,) =x» + tt£(x m ) = x» + A%(x,) + (4) 

where 7r^(a; M ) behaves in the two-site model denned by GCi and GCj as a Nambu-Goldstone 
boson, which is eaten to provide the longitudinal component of a massive graviton [11]. In 
Eq. (jl]), the Nambu-Goldstone bosons 7r^(a; M ) have been rewritten in terms of vector fields 
and scalars <pji, where the (ftji provide the longitudinal components of the massive gravitons. 
In fact, transforming to the stationary (or unitary) gauge Y-\{x^) = x M , we observe that 
the "hopping" terms in Eq. ([3]) produce graviton mass terms of the Fierz-Pauli type [25]. 
By analogy with the Eguchi-Kawai model, the model with two discrete gravitational extra 
dimensions will therefore lead to a multi-graviton theory, containing one zero-mode graviton 
and a phonon-like spectrum of massive gravitons with masses of order ~ l/(Na), which 
simulates in the IR a linear tower of KK excitations. 

Now, let us consider small fluctuations of the metrics g % = r/^ + around Minkowski 
space and express the gravitons as h % = e l2m ' n ^ N h n , where we have introduced the index 
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n = (m, no) with m, rz.2 = 1, 2, . . . and where we have suppressed for simplicity the Lorentz 
indices. Similarly, we write the scalar modes of the Nambu-Goldstone bosons as cj) 1 ^ = 
^^ e i2n-n/jv^^ where k — 1 (k — 2) if the vector i — j points in the y\ (1/2) direction. 
Inserting Eq. (J3|) into Eq. ([3]) and transforming to momentum space, we observe that the 
action contains a part, which takes in the limit tiCAfa form of the type 



S = d 4 x N 2 M 2 



I! 



fc=l,2 ^ ' 



+ -.., (5) 



where we have omitted terms involving the transverse spin-1 modes and higher-order terms 
in the spin-2 excitations. After going to canonical normalization, h n — > h' n = NMh n and 
— fin = ("i+»2)M 0fc^ we g nc j f rom term in Eq. (jSJ), when evaluated for the 

lowest lying modes, the strong coupling scale 

A = {M Pl Nl&yi\ (6) 

which is similar to the corresponding value in a single discrete extra dimension. Now, in 
order to avoid effects that are non-local in theory space [12,13], the inverse lattice spacing 
must always be smaller than the maximal short-distance cutoff A max ~ a/ M§d / R. Assuming 
a compactification scale of the order of 1/R ~ 10~ 2 eV ~ (10 /zm) -1 and Mqd ~ 10 TeV, 
we thus obtain a cutoff A max ~ 0.3 MeV, which is by many orders of magnitude smaller 
than Mqe>. From a field theory point of view, it would therefore be necessary to provide an 
understanding of the hierarchy a" 1 < A max M 6D , when the lattice spacing a is dynamically 
generated in a possible UV completion of the theory of massive gravitons. 

This hierarchy will be even more enhanced in strong gravitational fields, where the higher- 
order terms omitted in Eq. (jSJ) can no longer be neglected. Generally, for a macroscopic 
body with mass M', the cutoff gets modified as A max -> A' max ~ A max (M Pi /M / ) 1 / 3 [11]. 
Currently, the most precise Cavendish-like tests of Newton's law use test masses of the order 
M' ~ (1 — 10) g (see, e.g., Ref. [26]). If we are only interested in local theory space models of 
gravity that admit meaningful predictions in such laboratory experiments, then the inverse 
lattice spacing must be even smaller than a cutoff A max ~ 1 keV. This suggests to consider 
theory space formulations of large extra dimensions which naturally generate small inverse 
lattice spacings. Going to the extreme limit where the inverse lattice spacing becomes as 
small as a few eV, such models would have the additional benefit that the structure of theory 
space becomes potentially accessible in low-energy experiments via neutrino oscillations, 
when a right-handed neutrino is propagating in the latticized bulk. 

Note that, to prevent the SM from getting strongly coupled at ~ A max , we would require 
a perturbative description of the gravitons above this scale. In absence of such an UV 
completion for gravity, however, we shall in the next section analyze instead a deconstructed 
£7(1) gauge theory, which provides an UV completion for the theory of massive gauge bosons. 
This has the advantage that one can readily formulate an explicit mechanism which produces 
sub-mm lattice spacings from energy scales above the lattice cutoff, while ensuring that the 
dynamics of the SM fermions always remains perturbatively sensible. 
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Figure 2: Moose diagram for the deconstructed U(l) gauge theory on a two-dimensional disk. 
Each circle corresponds to one U(l)i = Gi (i = 0, 1, 2, . . . ,N) gauge group. An arrow pointing 
toward (outwards) a circle denotes a field with negative (positive) charge under this group. The 
link fields Qi t i+i define the boundary of the disk, while the radial links Qo,* connect the gauge 
group in the center with the sites on the boundary. 

3 Deconstructed U(l) on a disk 

In this section, we will study the deconstruction of a U(l) gauge theory in a non-gravitational 
extra dimension compactified on the boundary of a two-dimensional disk. This theory space 
has been analyzed in the context of supersymmetry breaking [7] and the doublet-triplet 
splitting problem [27]. Various properties of supersymmetric deconstructed U(l) models 
have also been studied in Ref. [28]. On the boundary of the disk, sub-mm lattice spacings 
can be generated from a hierarchy between the masses of the link fields. The resulting coarse 
latticization is experienced by a right-handed neutrino, which propagates on the boundary 
of the disk and mixes with the SM neutrinos located in the center. The deconstructed U(l) 
acts on the SM fermions as a gauged B — L symmetry that is broken at the TeV scale. 

3.1 Gauge sector 

Consider the deconstruction of a U(l) gauge theory, which is defined on the boundary of 
a two-dimensional disk. Our deconstructed theory is described in four dimensions by a 
jj^N+i = U^L U(l)i product gauge group, where each gauge group U(l)i corresponds to a 
site in theory space. The N+l sites are connected by 2N scalar link variables Qo,i and Qi^+i 
(i = 1, 2, ... , N), each of which carries under exactly two neighboring gauge groups U(l)i 
and U(l)j the U(l)i x U(l)j charges (+1, —1) and transforms trivially under all the other 
gauge groups. Specifically, for i — 1, 2, . . . , N — 1, a link Q^i+i is charged as (+1, —1) under 
the product group Z7(l)j x U(l) i+ i, while the link Qn,i = Qn,n+i carries the U(1)n x U(1)i 
charges (+1, —1). Thus, the sub-graph defined by the link fields Qij+i has the geometry of a 
latticized circle. Each link Q Q i (i — 1, 2, . . . , N) carries the Z7(l) x U(l)i charges (+1, — 1), 
which leads to a "non-local" theory space, since any two sites are connected through at most 
two links. The theory space of this model is conveniently represented by the "moose" [29] 
or "quiver" [30] diagram in Fig. [2l The gauge group Z7 (l)o corresponds to the center of 
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the disk, while the other gauge groups £7(1)^0 define the sites on the boundary and are 
connected by the "boundary links" Qii+i. The "radial links" Qo,i connect the center to the 
sites on the boundary. We wish to reiterate that, in contrast to the previous section, 4D 
gravity is simply added to the 4D model, i.e., our deconstructed extra-dimensional manifold 
is non-gravitational. 

It is useful to consider the global symmetry which corresponds to a 2tt/N rotation of the 
disk and acts on the link fields by 



Q 



N,l 



Qi,2j Q 



Q,N 



Qo,l: Qi,i+1 ~ * Q 



i+l,i+2) 



Qo,i ~^ Q 



0,i+l) 



(7) 



where % = 1,2, . . . , N — 1. Using the gauge degrees of freedom, we can always establish an 
equivalence relation Q12 ~ Qi,i+i, for i = 1,2, ...,N, which "identifies" the links on the 
boundary. Now, for the lattice gauge field "living" on the disk, the holonomy around each 
plaquette in Fig. [2] is trivial in the lowest energy state. As a consequence, the Wilson lines 
will break the symmetry U(1) N+1 x Z N down to f7(l)di ag x F, where C/(l)diag is the diagonal 
subgroup of U(1) N+1 and F can be taken as a diagonal product (linear combination) of a 
[7(1)0 g au g e transformation and the global Z N symmetry [27]. 

Initially, the scalar sector possesses a global U(1) 2N symmetry, which is broken by the 
gauge couplings, such that only a U(1) N+1 subgroup is preserved as the gauge symmetry 
of the model. Wilson line breaking of the global symmetry leads to 2N Nambu-Goldstone 
boson fields. At the same time, the Wilson lines break also N generators of the U(l) N+l 
gauge symmetry, which produces N massive spin-1 vector states by eating N of the Nambu- 
Goldstone boson fields via the Higgs mechanism. Thus, we are left with N classically massless 
Nambu-Goldstone bosons in the low-energy theory, which can, however, acquire a mass at 
the quantum level, since the large global U(1) 2N symmetry is explicitly violated in the gauge 
sector. 

The masses of the gauge bosons receive a contribution generated via the Higgs mechanism 
from the kinetic terms of the link fields ~ J2iLi(\F > i J ,Qo,i\ 2 + \F>^,Qi t i+i\ 2 ), where D^Qji = 
(<9 M + igjAjp — igiAi^Qj^ for e {(0, i), {%,% + 1)} is the covariant derivative, in which 
gfj and A ifM denote the gauge coupling and the gauge boson of the group Z7 (1)*- In the basis 
(Aq, Ai, A%, . . . , Ax), the (N + 1) x (N + 1) gauge boson mass squared matrix M\, which 
results from these terms after spontaneous symmetry breaking (SSB), is therefore 



M\ 



2 2 
9 v 



( N 
-1 

-1 
-1 

V : 



1 






1 






1 



1 2 2 

+ g u 







/o 

2-1 
0-12 

0-1 

V ; : 








where we have assumed, for simplicity, universal gauge couplings gi = g, while v and u 
denote the universal VEV's v = (Qo,i) an d u = (Qi t i+i), of the radial and the boundary link 
fields, respectively. After diagonalization of the mass squared matrix in Eq. (JSJ), we arrive 
at the gauge boson mass spectrum 



Ml 



0, Ml = g 2 v 2 + 4#Vsin 



,2„ 2 



, 7m 
It'' 



M 2 N 



2„.2 



(N+l)g 2 v 



(9) 
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where n = 1, 2, . . . , N — 1. We observe that this spectrum contains a zero mode, which 
would correspond to an unbroken {7(l)di ag (this symmetry will be broken later, when we 
introduce fermions). In Eq. (Q, the tower of mass squares M 2 , where n = 1,2, . . . ,N — 1, 
reproduces for n <C iV a spectrum of KK modes of the order (n/R) 2 = (ngu/N) 2 , which 
has been shifted to higher values by an additional universal contribution of the order (gv) 2 
provided by the radial links. One gauge boson with mass decouples at low energies 
for N — > oo. Thus, in this limit, the theory of massive gauge bosons becomes an effective 
description of a latticized flat fifth dimension with an inverse lattice spacing u (generated by 
the boundary links) and one lattice scalar (represented by the radial links), which acquires 
a VEV v in the 5D bulk. In the process, the link fields break the total U(l) product gauge 
group down to the diagonal subgroup U(1) N+1 — > £/(l)dia g , which can be further broken by 
suitable scalar site variables in the center that acquire a VEV at the TeV scale and will thus 
correspondingly modify the gauge boson spectrum in Eq. (jUJ). 

3.2 Large lattice spacings 

In order to determine the actual vacuum structure in more detail, let us now consider the 
scalar potential of the link fields in isolation. The most general gauge invariant renormaliz- 
able scalar potential of these fields then reads 

N 

v = X) [ m 'l^| 2 + M 2 \Q iyi+l \ 2 + fJ,Qo,iQi,i+iQl i+1 + A**Qo,i+iQi li+ iQo,i 
i=i 

1 1 i- r 

+ 2^ 1 l < 5°> i l 4 + 2^ 2 l^ i . i + 1 l 4 + ^s\Qhi+i\ 2 2^1 l^o.il 2 + ^4 l^5o,i| 2 X^ l^ojl 2 

+^s\Qi,i+l\ 2 X^ l^i.i+il 2 + (^6 < 5o,j < 5i,j+l < 5i+l,i+2<5o ) j + 2 + h.C.)j , (10) 

where the parameters m,M, and \i have mass dimension +1, while Ai, A 2 , A^ 7 , A4 , and \§ 
are dimensionless real parameters of order unity and A6 is a complex-valued order unity 
coefficient. Note that the potential V is invariant under the global symmetry in Eq. (J7J). 
From the point of view of usual effective field theories, the dimensionful parameters m, M, 
and fi may take any value in the UV desert between ~ ITeV and Mpi. We will consider here 
the interesting case where these masses exhibit a hierarchy m fi ~ M, i.e., the boundary 
link fields are much heavier than the radial link fields. To be specific, we assume that the 
mass m is close to the TeV scale, i.e., m ~ 1 TeV, while [i and M are of the order the 
usual B — L breaking scale \i ~ M ~ Mb-l — 10 15 GeV. Note that an understanding of 
the smallness of the parameter m with respect to Mb-l m &y require mechanisms similar to 
those who give a solution to the //-term problem in supersymmetric theories and will not be 
specifically discussed here. 

To explicitly minimize the scalar potential V , we shall now make the simplifying assump- 
tion that the parameters /1 and A6 in Eq. (fTOj) are real. Actually, in a supersymmetric case, 
holomorphy of the superpotential sets A 6 — ► and one could rotate the phase of // into the 
Yukawa couplings of the neutrinos. Therefore, we argue that our basic results concerning the 
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magnitude of the VEV's will not be significantly altered, when considering the more general 
case of complex \i and A 6 . Now, taking m 2 < and fi < 0, while M 2 > 0, the potential 
V in Eq. fflUj) has an extremum [35], which is given by u = (Q i)i+1 ) and v = (Qo,i), for 
i = 1, 2, . . . , N, where u and v are real and equal to 

m s 

U ~ 2[\ 1 + {N-l)\ ^ )MZ-^ , 1 } 

z.e., the boundary links and the radial links respectively acquire universal VEV's. From 
Eqs. (11 II) we observe that u and v 2 become ~ 1/N suppressed in the large N limit. However, 
let us now consider the opposite situation, where N ~ O(10), i.e., the number of sites is 
kept moderate or small. In this case, we observe that the choice of mass scales m ~ 1 TeV 
and /i ~ M ~ Mb-l — 10 15 GeV generates for the boundary link fields Qi^+i a small VEV 
of the order u ~ 10 _1 eV, while the radial link fields Qo,? acquire an unsuppressed TeV scale 
VEV v ~ 1 TeV. In other words, the model generates from mass scales in the UV desert 
of conventional 4D theories an inverse lattice spacing u ~ (/im) -1 in the IR desert of large 
extra dimensions. The suppression of u due to the hierarchy m < /i ~ M is similar to the 
type-II seesaw mechanism [36] and, in fact, the structure of V in Eq. ( TTUT) can essentially 
be viewed as a replication of the model in Ref. [37]. It is the replication of gauge groups 
on the boundary, which allows here to interpret ~ u^ 1 as the sub-mm lattice spacing of a 
deconstructed large extra dimension. 

Note that our mechanism for the generation of sub-mm lattice spacings differs from the 
model in Ref. [35] essentially in the choice of the dimensionful parameters m, fi, and M. In 
Ref. [35], they are of the orders m ~ /x ~ 10 2 GeV and M ~ 10 9 GeV. Having /j ~ M, 
however, would be an automatic consequence in a supersymmetric version, where the tri- 
linear plaquette terms in Eq. ffTOl can emerge from the F-terms of the superpotential. 



3.3 Inclusion of fermions 

In our U(1) N+1 model, we will first extend the three generations of SM fermions by three 
fermions Ni, N 2 , and A^, which are singlets under the SM gauge group G$m- Then, the 
three fermion generations are put on the center of the disk by assuming that they carry 
nonzero U(1)q charges, but are singlets under the other gauge groups [7(1)^0- Note that 
the addition of matter fields as site variables on the center leaves the Zn symmetry in 
Eq. (J7|) unbroken. We suppose that the leptons £ a = (v a , e a ) T and e c a (a = 1, 2, 3 is the 
generation index) are charged under U(l)o as +1 and —1, respectively, while the quark 
doublets q a = (u a , d a ) carry a Z7 (l)o charge —1/3, and the isosinglets u c a , and d c a are given 
the U(1)q charges +1/3. The SM singlets N\, N 2 , and N 3 , carry the U(l) charges —4, —4, 
and +5, respectively. Since the U(1)q charges of the SM quarks and leptons are identical 
with their B — L quantum numbers, it is easily seen that the model will be free from axial- 
vector [31] and gauge-gravitational [32] anomalies. Notice that this is slightly different from 
the usual way of gauging B — L, where three right-handed SM singlet neutrinos carry a B — L 
charge —1 [33]. With our charge assignment, however, the Yukawa couplings of the active 
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neutrinos to the fields N a are suppressed by many powers of Mpi and are thus negligible. 
Suitable SM singlet scalar fields S with nonzero U(l) charges can then allow renormalizable 
Yukawa couplings ~ SN a Np and break the U(l)o symmetry around the TeV scale. (For a 
recent detailed analysis of breaking B — L at the TeV scale see, e.g., Ref. [34].) The fields 
N a , which were only introduced for the purpose of anomaly cancellation, will then decouple 
below a TeV. 

Next, we include a SM singlet fermion, which appears with respect to the SM interactions 
as a right-handed neutrino propagating on the boundary of the disk. In the deconstructed 
space, the bulk fermion is represented by N right-handed neutrinos (i = 1,2, . . . , N), 
which are put as site variables on the boundary of the disk. Here, has a charge —1 under 
the group U(l)i, but is a singlet under the other gauge groups [/(l)^. In the Weyl basis, 
we can decompose each field as ^ = (vm, v Ri ) T , where i/m and u Ri are two-component 
Weyl spinors. The mass terms of the neutrinos will be discussed in the next section. 



4 Neutrino masses 

In this section, we will analyze the kinetic term of the latticized right-handed neutrino and 
its mixing with the SM neutrinos. Unwanted higher-dimension operators can be eliminated 
by refining the triangulation of the disk. 



4.1 Latticized right-handed neutrino 

The deconstructed model presented so far describes a non-local theory space where any 
two sites are connected by at most two links. This is in contrast to the continuum theory 
for neutrinos in large extra dimensions [17,18], where a local interaction of a massless right- 
handed neutrino in the bulk leads to a ~ [Mp\)~ x suppressed coupling to the active neutrinos. 
Moreover, the model is in its present form vector-like and allows unprotected Dirac masses 
~ Mp\VpiV c Ri . However, if we treat the latticized right-handed neutrino as a massless Wilson 
fermion [38] propagating on the boundary of the disk, which has lattice spacings in the sub- 
mm range, we can have only small Dirac masses ~ uvpiV c Ri) where u ~ (/im)" 1 is the inverse 
lattice spacing. In order to remedy this problem and make contact with the 5D continuum 
theory in Refs. [17,18], we introduce for each site on the boundary of the disk a pair of scalars 
Xn and (f) n (n — 1, 2, . . . , N), which are Gsm x ^(l)^ 1 singlets, and assume a discrete Zqm 
symmetry (M is an appropriate integer) acting on the fields as 



v Rn - 

Xn 
Qo,n 

e c 

K 

N h2 



-i27r(2n+3)/M 
c An? 

- ^ /M Qo,n, 

ja-K/M c 

e -i27r/(3M) c 



48tt 



/M Ni 



2- 



V 



Rn 



-i27r(n+l) 2 /M 



Rni 



. e i2?r(2n+3)/(2M)0 

-i2-K/Mf 

i27r/(3M) 
c Hai 
e -i27r/(3M)^c 

>e -iWn/M N ^ 



(12) 



where n — 1, 2, . . . , N and a — 1, 2, 3 runs over all three generations. Note that the left- and 
right-handed SM fermions carry opposite charges under the Zqm symmetry, and hence, the 
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Yukawa couplings of the quarks and charged leptons will remain unsuppressed. Furthermore, 
we note in Eq. ffl2l that the potential V in Eq. fflUj) remains invariant under the Zqm 
symmetry. In Appendix we show that by adding extra fermions on the boundary the 
Zq M symmetry can be promoted to a discrete gauge symmetry, which would be protected 
from quantum gravity corrections [39,40]. It turns out that in the effective theory (ignoring 
the enlarged gauge symmetry at high energies) all dangerous triangle diagrams would add up 
to zero. It is interesting to note here, that the SM possesses an anomaly-free Zq symmetry 
which can ensure nucleon stability for new physics scales as low as ~ 10 2 GeV [41]. 

In the deconstructed theory, the action for neutrino masses which includes all renor- 
malizable interactions with 0(1) Yukawa couplings and the most general U(l) N+l invariant 
dimension-five operators consistent with the discrete Zqm symmetry, can now be written as 

^mass ^wilson ^mt ^dim5; (l*^) 

in which the different parts are given by 

f N Q X 

Wilson = / d 4 X ^ UV Rn ( n ^ +1 *4 (n+1) ~ ^*4n) + h " C - » ( 14a ) 
n=l 

St» = J d*x^£ a eHQ* 0>1 v m + h.c, (14b) 
f N Y n 

SdimS = d 4 xYjfQl,nQo,n+lVRnV R{n+1 )+h.C., (14 C ) 
w 1 J 



n=l 



where e = icr 2 contracts the SU(2) indices, while Y a and Y n are (complex) dimensionless 
(9(1) Yukawa couplings. The non-renormalizable operators <S m ^ and <Sdim5 are generated 
at the string or "fundamental" scale Mf ~ (10 17 — 10 18 ) GeV, where a value as low as 
Mf ~ 10 17 GeV could be understood in M-theory [42]. In Eq. (I14al) . let us assume that 
the Xn acquire a universal VEV (x n ) = u {n = 1, 2, . . . , N), which is equal to the inverse 
lattice spacing u defined by the universal VEV's of the boundary links in Eq. filial) . We 
will comment on a possible origin of the order of this mass scale for the VEV's of the Xn 
later on. With the identification (x n ) = u, the action 5 wi i son in Eq. f)14ap takes the form 
of a Wilson-modified latticized 5D kinetic term that describes the propagation of the right- 
handed neutrino on the boundary of the disk, which is interpreted as a fifth dimension. The 
fields u RQ and v R ( N+1 ^ in Eq. (I14ap are determined by u c Rl and z/f^ only up to a discrete Z 2 
"gauge transformation" reflecting the topology of the disk. The Z 2 symmetry is associated 
with the existence of non-trivial or twisted field configurations [43] for the latticized right- 
handed neutrino, which are characterized by distinct spectra in the low-energy theory. In 
Eq. f!14ap. we define (v R0 , v c m ) = ±(v RN , v L RN ) and (v R (n+i), *4(tv+i)) = M u Rh *4i)> where 
"±" distinguishes between twisted (—1) and untwisted (+1) fields. The effects of twisted 
field configurations in deconstruction have been extensively discussed in Ref. [44]. 

Upon using the mechanism in Sec. [3] for generating small inverse lattice spacings u ~ 
(/im) -1 , the latticized 5D kinetic term *S w iison leads then to an effective action for KK modes 

f N 

Skk = / d*x Y uv Rn (v c R{n+l) - z4„) + h.c. (15) 

n=l 
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In the basis spanned by (vri, ur2, ■ • • , vrn) and 
defines a Dirac mass squared matrix M 2 , which explicitly reads 



V R\i V R2i 



RN 



) , the action iS>kk in Eq. f|T5|) 



M 



it 



V 



-T 



-1 
2 



-T \ 



where T = ±1 and the blank entries are all zero, 
are found to be the eigenvalues of the matrix M 2 
and untwisted (T = +1) fields the mass spectra 



The squared masses m 2 of the fermions 
Thus, we obtain for twisted (T = — 1) 



4-u 2 sin 2 



n 



1/2)tt 



(twisted), 



4u 2 sin 2 



n 



1)tt 



V 



(untwisted), (16) 



N 

where u ~ (/im) -1 and n = 1, 2, . . . , N. We hence observe that Skk reproduces in the IR 
for VR n and v Rn always a tower of KK excitations with ~ (/mi) -1 masses, which becomes 
for large A" indistinguishable from the lightest KK modes of a right-handed bulk neutrino 
in sub- mm sized continuum extra dimensions. Note in Eq. (TIB"]) , that a zero mode is absent 
for twisted fields. 

In the above discussion, we require that the Xn acquire the small VEV (x n ) = u ~ (a™) -1 
to allow the interpretation of 5 w iison as the Wilson action for a massless right-handed neutrino. 
This energy scale has been generated for the VEV's of the boundary links by the mechanism 
in Sec. |3j Since the potential for the Xn and (p n is qualitatively similar to the potential 
V of the link fields in Eq. ffTUl) . a variation of this mechanism can also produce the right 
energy scale for when Xn and <p n take the roles of the boundary and radial links, 

respectively. For this purpose, we suppose that the (f> n have masses around the TeV scale 



have masses M x of the order the Planck scale M x 



M 



/• 



~ 10 2 TeV, whereas the x 
Additionally, we take in the renormalizable ZgM-invariant interactions ~ f^n^n^nXn the 
dimensionful couplings jl n to be jl n ~ Mf. By the same arguments as in Sec. [3l we find that 
the Xn can acquire a VEV in the range (x n ) — m 
inverse lattice spacing u in Eq. filial) . 



\/Mf ~ (jum) 1 , which is of the order the 



4.2 Non-renormalizable terms 

The interaction of the active neutrinos v a with the right-handed neutrinos on the boundary 
of the disk is introduced by Sf^ in Eq. (]14bj) . In theory space, we identify the dimension- 
five term [45, 46] in Sf^ with a Wilson line type effective operator, which connects the 
active neutrinos (in the center of the disk) with Vrx (o n the boundary) via the link Qo,i 
(see Fig. [3]). Let us now go to the basis where iSkk i n Eq. f[T5"j) is on diagonal form and 
consider the lowest lying mass eigenstate i/ R belonging to Skk- After setting Qo,i to its VEV 
(Qo,i) — v — 1 TeV, the Wilson line type operator generates an effective Yukawa interaction 
~ Y a £ a eHiy' R v/(\/ r NMf), which is suppressed by a factor ~ v/{y~NMf) with respect to the 
electroweak scale. For a string scale Mf ~ 10 17 GeV and small A^ we thus obtain Dirac mass 
terms m£) a is a is R , with Dirac masses rriDa = Y a (H)v/(\fNMf) ~ 1CT 3 eV. 
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<Qo,i) 




(Qo,n) 



VR( n +\y 



(Qa,n+\) 



(a) (b) 

Figure 3: Non-renormalizable dimension- five operators generated at the fundamental scale Mf. 
When the link field Qo,i acquires a VEV (Qo,l) — ITeV, this operator generates for Mf ~ 10 17 GeV 
a Dirac mass of the order 10 -3 eV, which mixes the active neutrinos v a with the tower of KK states 
via an interaction with Vr\. In theory space, this dimension-five term corresponds to a Wilson 
line type effective operator, which connects the SM neutrinos v a in the center with ur\ on the 
boundary of the disk via the link field <5o,i- The operators in (b) generate Dirac masses of the 
order 10 _2 eV for the right-handed neutrinos and are neglected with respect to the nearest neighbor 
hopping terms of the order 10 -1 eV. 



It is instructive to compare the effective Yukawa interaction generated by with the 
5D ADD scenario. Here, a right-handed bulk neutrino vr couples to the active neutrinos on 
the SM brane through a local interaction [17, 18] 

= J d 4 x^e a (x)eH(x)Mx,y = 0), (17) 

where y is the coordinate along the fifth dimension compactified on a circle with circumfer- 
ence 2%R, the coefficients are dimensionless 0(1) Yukawa couplings, and the SM lepton 
doublets £ a as well as the Higgs doublet H are 4D fields trapped at y = on the SM brane. 
Note that while vr\ in Sf£ has mass dimension 3/2, the 5D fermion vr in Sf® has mass 
dimension 2. After expanding v R in KK modes as v R (x, y = 0) = (2txR)~ 1 ^ 2 J2 n u Rn( x ) an d 
using the relation 2nR = M Pl /M^, it is seen that the interaction in Eq. (ITT)) gives rise to a 
Dirac type coupling ~ Y® £ a 6Hi>RQM* / Mpi between the active neutrinos and the zero mode 
uro, which is M^/Mpi suppressed. Since v ~ M* ~ 1 TeV, we thus find that, in the limit 
of coarse latticizations y/N ~ (9(1 — 10), the couplings between the active and the right- 
handed neutrinos generated by and become suppressed by factors of similar orders 
v/{yfNMf) ~ M m /Mp\. However, despite this numerical coincidence, the two models differ 
in an interesting way: while the smallness of the Dirac type coupling in Sf^ emerges from a 
volume suppression factor (i.e., from the large number of KK modes below M*), the small 
Dirac mass generated by is rather a result of the separation between the site where the 
SM fermions are located and the boundary of the disk as compared to the length scale Mj 1 . 

The dimension-five operators contained in iSdi m 5 hi Eq. fl!4cj) give for Mf ~ 10 17 GeV 
rise to Dirac mass terms between VR n and ^( n+1 ) that are of the order 1CT 2 eV (see Fig. [3]). 
However, since vr\ is the only right-handed neutrino which couples "directly" (at the non- 
renormalizable level) to the active neutrinos, we may treat the terms in iSdims as subleading 
corrections to <S w ii son , which gives Dirac masses of the order 10 _1 eV, and ignore them in 
the following discussion. Let us now, instead, consider a more attractive possibility to 
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Figure 4: Spider web theory space for the suppression of unwanted higher-dimension operators. 
The inner disk is defined like in Fig. [2] and the latticized right-handed neutrino propagates on the 
outer circle formed by N additional U(l) gauge groups J7(l)^ = G\ (i = 1,2, . . . , N). The SM 
fermions are placed on the site corresponding to G% (gray site). The link fields Q\ i+l on the outer 
circle have a mass of the order 10 12 GeV, whereas all other link fields have masses of the order 
10 2 TeV. Due to the plaquette-structure near the boundary involving four links, the operators 
analogous to Fig. [3] (b) are suppressed by an extra factor ~ 10~ 12 . 



suppress the unwanted non-renormalizable terms of the type shown in Fig. E^b) by making 
only use of the plaquette-structure of the model. For this purpose, we will assume that the 
non-local theory space introduced in Sec. [3] is actually part of a larger "spider web theory 
space" [7] as shown in Fig. HI This theory space is obtained from the disk in Fig. [2] by 
adding N extra U(l) gauge groups (i = 1,2, . . . , N), which gives the total gauge 

group U(1) 2N+1 = Uf =0 U(l)i x IFf =1 [/(l);.. Each pair of factors £/(l)J and U{l)i is connected 
by a link field Q' ii which is charged under U(l)i x £/(l)^ as (+1, —1) and is a singlet under 
the other gauge groups. Two neighboring groups £/(l)^ and U(l)' i+1 , where i ~ i + N, are 
connected by a boundary link field Q' ii+1 , that is charged as (+1, —1) under C/(l)^ x U(l)' i+1 
and transforms trivially under the other gauge groups. 

In analogy with Sec. [31 we suppose that the latticized right-handed neutrino propagates 
in Fig. H]on the outer circle defined by the links Q' ii+ i- Contrary to the previous non-local 
theory space example, however, the SM fermions are now placed on the site associated with 
G\ on the inner circle (see Fig. [4]). The £7(1) i charge assignment of the SM fermions is 
similar to that in Sec. [3] with U(l) replaced by Z7 (l)i- Next, we suppose that the boundary 
link fields Q[ i+1 on the outer circle have a common mass of the order of an intermediate 
scale 10 12 GeV and positive mass squares. The remaining link fields Qo t i, Qi,i+i, and Q' i{ , 
on the other hand, are supposed to have masses of the order 10 2 TeV and negative mass 
squares. By the same arguments as in Sec. [3], we then find that the corresponding scalar 
potential is extremized for VEV's (Q' ii+ i) — 1 eV, while all other link fields have VEV's of 
the order (Qo,i) — (Qi,i+i) — (Q'a) — 10 2 TeV. The separations of the sites on the outer 
circle are therefore in the sub-mm range, whereas the inverse lattice spacings between all 
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other neighboring sites are of the order 10 2 TeV. As a result, the active neutrinos couple to 
the latticized neutrino on the outer circle via an analog of the operator in Eq. f!14bj) . 
where only Q ,i has been replaced by Q'n- More importantly, the "dangerous" dimension- 
five operators of the type shown in Fig. [3] (b) are now replaced by dimension-six operators 
as v Rn v c Rin+1) Q*Q )n Qo,n+i/M f -> VRnV R{n+1) Q'l n Qn,n+iQ' n+ i, n+ i/M 2 f . As a consequence of 
the plaquette-structure near the boundary of the spider web theory space, the unwanted 
higher- dimension operators are suppressed by an extra factor ~ 10~ 12 and become therefore 
irrelevant. We thus see, that the model in Sec. I4.1l is completely reproduced for a fundamental 
scale Mf ~ Mpi ~ 10 18 GeV, but now without the unwanted higher- dimension terms of the 
sort given in Eq. fll4cl) . 

Up to now, we have been considering the generation of Dirac neutrino masses in a model, 
where B — L is preserved by the link fields. In order to understand recent solar [47], at- 
mospheric [48], reactor [49], and accelerator [50] neutrino data we assume - contrary to the 
usual type-I seesaw mechanism [51] - that B — L is broken at the TeV scale. One possibility 
is here provided by versions of the Babu-Zee model [52] , which can be easily implemented 
in our model to generate radiatively Majorana neutrino masses locally on the site where 
the SM fermions are located. However, in what follows, we will not further specify the de- 
tailed mechanism which generates the Majorana masses of the usual neutrinos. Instead, we 
will always assume the presence of suitable Majorana mass terms and concentrate in the 
deconstructed U(l) model on the mixing of the SM neutrinos with the KK modes that is 
introduced by the Dirac neutrino masses. 

5 Mixing with Kaluza— Klein modes 

In this section, we will consider the neutrino mass and mixing terms of our model by spe- 
cializing to the simplifying case of only one single active neutrino v coupling to the latticized 
right-handed neutrino, which we treat as a Wilson fermion. From the action in Eq. ([13]) . 
we thus obtain in this case after SSB for the relevant neutrino mass and mixing terms the 
action density 

N-l 

C v m = m v vv + \[~Nm D vv m + uv RN {Tv c Rl - v c RN ) + VR n {v c R{n+l) - v c Rn ) + h.c, (18) 

n=l 

where the small inverse lattice spacing u ~ (/mi)" 1 has been generated by the mechanism 
in Sec. |3j and the parameter T = ±1 describes a twisted/untwisted right-handed neutrino. 
In Eq. f fl8|) . the Dirac mass type coupling y/Nrrif) = Mj 1 (H)(Q ^) ~ 10~ 2 eV arises from 
the higher-dimension operator shown in Fig. [3J (a) (or the analogous term in the extension 
to spider web theory space with Q 0>1 replaced by Q[ 1; see Sec. @J, while the Majorana mass 
m u ~ 10 -2 eV has some other origin and may, e.g., emerge from a radiative mechanism as 
mentioned in Sec. HI For convenience, we have chosen for the second term in Eq. (TlBl a nor- 
malization factor x/iV, which is related to the volume suppression factor in the corresponding 
5D continuum theory. 

The action density £^ in Eq. ffl8l) translates into a (2N + 1) x (2N + 1) neutrino mass 
matrix from which we determine by diagonalization the neutrino mass and mixing param- 
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eters for the different cases iV odd/even and T twisted/untwisted. Since m v , yNmn <C u, 
it is useful to define the quantity e = \/~Nm£>/u < 1 as an expansion parameter in per- 
turbation theory and diagonalize the matrix MM^ in several steps. First, we bring the 
latticized fermion kinetic term in Eq. (|18|) on (approximately) diagonal form by applying a 
transformation MM^ — > / M M I with a suitable unitary matrix U. The mixing matrices 
U for the different possible cases are explicitly given in Appendix [B] For definiteness, let us 
consider the case T = — 1 and N even, the other cases follow then in similar ways. Trans- 
forming to momentum space with respect to the latticized dimension defines a new basis 



[V, VRN/2+U 



vrn, Vri, v R2 , 



V 



RNSi 



in which the resulting mass squared matrix reads 



M 



u 



V 



A 

7 

7 
7 
a x 

a-N/2 

^N/2+l 
1>N 



^N/2+l + <5 



7 

6 



7 
6 



a N/2 ON/2+1 



>N 



\ 



Xn~i +3 5 
S \ N + 5 



Ai 



A 



JV/2+1 



Xn j 



where the blank entries in this matrix are all zero. The nonzero elements are given by 



7 



6 



2m 2 D 

9 : 

U 



VzmD (2n — 1W 
sm 



u 



N 



(19) 



(20) 



for n = 1, 2, . . . , N/2 and 



b„ 



u 



•1 + cos 



(2n- 1)tt 



(21) 



for n = A^/2 + 1, N/2 + 2, 

A = 



. . . , N. In Eq. (Tl9l . the masses A and and A n are 

Nm 2 D + ml . 2 (n - 1/2)tt 

— and A„ = 4 sm - 



iV 



(22) 



for n = 1, 2, . . . , iV. The coefficients A n show the characteristic doubling of KK modes of a 
phonon-like spectrum, since they satisfy A n = Ajv-n+i for n = 1, 2, . . . , N/2. In Eq. ( fl9|) . note 
that the KK states Om, VB2-, ■ ■ ■ ■> vrn/2 exhibit no Yukawa interaction with the active neutrino 
v, and hence, decouple completely from the SM interactions. Next, we apply to the basis 
{y, vrn/2+i, ■ ■ ■ , vrn, Vri, v R2 , ■ ■ ■ i ^rn) a sequence of rotations by defining the orthogonal 
states 

VRn = S n^R n + °npR(N-n+V\ aI1 d V Rn = C n^R n ~ S n^R{N-n+l)i (23) 



for n 



1,2, 



R{N-n+l) 

N/2. In Eq. ( J231) . we have s n = cos (""j^ 2 ) 71 " anc l Cn = _ s j n ^ — Lri^. 



N ^ ^ n - ^ N -. Crudely, 
this corresponds to "rotating away" in Eq. ( fl9l) half of the interactions ~ which reduces 
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the degeneracy of the problem from three-fold to two-fold. In the new basis, the mass squared 
matrix is given by 



/A 7 ■ 

7 + $ ■ 



M 



u 



7 
7 
d\ 

\d>N/2 



5 
S 



7 7 dt d 2 ... d N/2 \ 

5 5 



X 2 + 5 5 
5 \! + 5 



Ai 



(24) 



Aat/2/ 



where 



dr. 



2V2m D 



in 



sin 



1/2)tt 



(25) 



u N 

for n — 1,2, ... , N/2. Here 7, 5, A, and A n are the same as in Eqs. (I20I) and (|22|) . In the con- 
tinuum limit, we expect for large N to recover some relevant characteristics of a continuous 
large extra dimension. In order to match onto the 5D continuum theory, we will compare 
our model with the one given in Ref. [19]. Since in this model Majorana masses are absent, 
we assume in Eq. f l24|) that m u — > 0, which implies that 7 — > 0. Furthermore, the matrix ele- 
ments 5 are small in comparison with the quantities A n and can therefore be neglected when 
calculating to lowest order. This means that the N/2 states v c RX , v c R2 , . . . , v RN / 2 spanning in 
Eq. Q33D the top-left N/2 x N/2 submatrix with entries A< + 5 (i = 1, 2, . . . , N/2) on the 
diagonal, decouple from v. Consequently, we end up with just one KK tower of N/2 states 
u Rn ' (i = 1,2,..., N/2), which span the last N/2 rows and columns of M 2 in Eq. ([24]). The 
remaining entries in M 2 become for n <C iV asymptotically equal to 



u 2 d r . 



\f2rriD 
R 



(n - 1/2), m 2 A 



it A r 



1 ( 
R 2K 



l/2f 



(26) 



where we have used the fact that u = N/(2ttR). We will match our model onto the 5D 
continuum theory by setting mo = y(H)M*/Mpi, where y is some 0(1) Yukawa coupling in 
the ADD scenario. With this identification, our model reproduces for the case N odd and T 
untwisted (see Table [[]) in the IR exactly the effective neutrino mass squared matrix of the 
5D continuum theory for neutrino oscillations in extra dimensions as discussed in Ref. [19]. 

Next, we will diagonalize M 2 by using two- fold degenerate Rayleigh-Schrodinger pertur- 
bation theory. We start by rewriting the mass squared matrix M 2 as M 2 = Mg + eMf , 
where Mq is a diagonal matrix, M 2 is the perturbation matrix, and e = \/NmD/u ^ 1 
is the small expansion parameter. In order for perturbation theory to be valid, we require 



that \((t$\eMl\(tF>) / (E K "> - E) u >)\ < 1 for % ± j, where |$" ; ) denotes the zeroth order 



t(0)\ 



,(0) 



,(0) 



eigenvector, the corresponding eigenvalue, and r and s are the degeneracy indices. This 
means that we will require that 
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(0)\ 
ir 



2V2> 



mou sin 



(n-l/2)7r 
N 



Nm 2 D + m 2 — 4m 2 sin 2 — ^ 2 ' >7r 



N 



< 1 and 



y/2: 



m D m u 



Nm 2 D 



mf 



4M 2 sin 2(liz^ 



< 1. (27) 
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T 


N 




<^n=l,2,...,n max -l 




An=l,2,...,n max 


+ 1 
+ 1 
-1 
-1 


even 
odd 
even 
odd 


N/2 + 1 
N/2 + 1/2 

N/2 
N/2 + 1/2 


2V2m D r,[ n ( n ~ 1 ) 7r 
u " N 

2V2m D Sin ("-D- 
u N 
2V2m D ■ (n-l/2)7T 

M " JV 

2V2m D • (n-l/2)7T 

M " iV 


2mr> 

2V2rn D (V-I)tt 
« ljlU 2V 

2V2m D (iV-l)Tr 
u 2V 
2mj3 


4 sin v N ' 
4sin 2(n - 1) " 
4sin 2(n -^ 2 ^ 
4sin 2(n ^ 



Table 1: Relevant entries in the neutrino mass matrices of the type shown in Eq. (|24p in the 
limit 7^0 for the different cases T twisted/untwisted and N odd/even. The first and second 
line reproduce for n <C N exactly the 5D continuum theory results of Ref. [19] and the third line 
corresponds to the example considered in the text. In all cases, the parameter A is as given in 
Eq. ([22]). 



From these relations we note that perturbation theory will not be valid for arbitrarily large 
N, since the denominator becomes singular at some point when Nm 2 D + m 2 ~ 4m 2 sin 2 (n — 
l/2)n/N. For the other cases, T = —1, N odd and T — 1, N odd/even, one obtains 
essentially the same constraints. Now, we apply perturbation theory to this problem and 
obtain the matrix that diagonalizes M 2 as a result. We denote this matrix by W^ k \ where 
k denotes the fcth order in perturbation theory. Thus, the mixing matrix V which relates 
the original basis to the mass eigenstate basis via MM^ — > V T MM^V* is given by 



V = XJP X P 2 ■ ■ ■ P N/2 W {k) 



(28) 



where P n are the rotation matrices associated with the state redefinitions in Eq. ( 1231 . and 
W^ k \ as stated above, is the matrix of eigenvectors of M calculated to some order k in 
perturbation theory. The first row is what is of interest to us, since it gives the relevant 
mixing angles of v with the bulk modes. It will be entirely determined by W^ k '. Thus, to 
lowest order in perturbation theory, we find 



V 



(\ eA 1 

* 



eA N / 2 eBi 



(29) 



y* • ■ ■ * / 

which is an (N + 1) x (TV + 1) orthogonal matrix, where 

(n-l/2)7r 



A r . 



m, 



sin 



2 (n-l/2)7r _ A 
N 4 



and B r . 



sin 



N 



sin 



2 (n-l/2)7r 
N 



(30) 



for n = 1,2, ... , N/2 and the elements denoted by * are not relevant in the following dis- 
cussion. Note that one can diagonalize MM^ in other ways then the one described above. 
For example, one could have applied four-fold degenerate perturbation theory directly to the 
matrix MM^ . One can show that this gives the same result for the final neutrino oscillation 
probabilities. However, reducing the degeneracy makes the problem much easier to handle. 
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6 Neutrino oscillations 



Global analyses have well established that the standard active three-flavor neutrino oscilla- 
tions with mass squared differences of the orders of magnitude Am^ — 8.1 ■ 1CT 5 eV 2 and 
Amfu ~ 2.2 ■ 1CT 3 eV 2 are in excellent agreement with neutrino oscillation data (see, e.g., 
Ref. [53]). However, the KK modes of the latticized right-handed neutrino could provide 
a sizable subdominant component in solar and atmospheric neutrino oscillations, and thus, 
lead to new anomalies, which are in reach of more precise ongoing or future neutrino os- 
cillation experiments. In this section, we will derive the corresponding neutrino oscillation 
formulas for our deconstructed model. These will only be valid in the regime where the mix- 
ing parameters satisfy the constraints in Eq. (j27|) . For phenomenologically allowed values of 
the physical parameters, this means that the formulas below will in general be valid for a 
low or moderate number of sites, N < 10 — 100. 

In order to describe the neutrino oscillations in our model, we can write the flavor eigen- 
states as a linear combination of the mass eigenstates using the mixing matrix in Eq. (I29p . 
We find that 

/ N/2 N/2 \ 

\ n=l n=l / 

Here \uf) denote a flavor eigenstate for some flavor / G {e,fi, r} and \u), \u n ), and \u n ) 
denote the mass eigenstates. We have also introduced a normalization constant K, which 
follows from the condition \(vf\vf)\ 2 = 1. Thus, we find from Eq. (131]) that 

N/2 

K 2 = l + e 2 J2 K + B D ■ (32) 

n=l 

Next, in the transition survival probability Pff = P(vf — > vf) = \ (vf\vf{t)) | 2 , the time- 
evolved state \vf(t)) is given by 



1 . (Nm 2 D +ml) 



N/2 N/2 



v f (t)) = -e- i as ( \u) + eY^Aj^\u n ) + eY^Bj^\v n ) \ , (33) 



n=l n=l 



where the phases (p n and the mass-squared eigenvalues vnr n equal [cf. Eq. (TIB]) ] 

(Nm 2 D + ml-ml)t . 2 . 2 . 2 (n-l/2)n , nA . 
<Pn = - — and m 2 =4u 2 sm 2 ^ ^M-, (34) 

in which E is the neutrino energy. Using Eqs. (|3T]) and (133]) gives for the case T = — 1 and 
N even 

N/2 2 

i + e 2 j2( A i + Biy^ ■ 



P _ 1 



n=l 



(35) 



For the other cases one finds the transition probability expressions in similar ways, first 
starting by applying the matrix U\ for the case one considers and then by using a set of 
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rotations similar to the ones given by Eq. fl23l) . Next, one applies perturbation theory and 
finds the mixing matrix from which the transition survival probability expressions follows. 
Thus, for the case T = — 1 and iV odd we have 



1 

K 4 



where 



Pff - — 



! + e 2 \J2(A 2 n + B 2 n )e^+f1e 



n=l 



(36) 



mr 



N 4- 



Nm? n +m1 



2 ■ 



Similarly, for T = 1 and iV even we find 



ff 



1 



N/2 



.Y-2 



n=2 



(37) 



where 



2 2 
7717,16 



AT (A^m^ + 



and /3 is the same as in Eq. ( 1361) . Note the presence of a zero mode in the phase fa. Finally, 
for T = 1 and A" odd we find 



Si 



1 



JV+l 
2 



n=2 



(38) 



where a is the same as in Eq. (13T1) . For the cases T = — 1, A n and S„ are given by Eq. (1311 
and the phases </> n are given by Eq. (13~4"|) . whereas for the cases T = 1 we have 



uV8N 



sm 



2 (n-l)Tr _ A 
N 4 



and B r . 



sm 



(n— l)7r 
TV 



sm 



2 (rc-l)7r _ A 
N 4 



(39) 



In this case, the phases are given by Eq. fl34|) . but with the masses [cf. Eq. (JTB 

(n — 1)tt 



2 a 2 ■ 2 

m„ = Au sm 



(40) 



In Figs. [5HU we have illustrated the different neutrino transition survival probabilities in 
vacuum as functions of L/E for the different cases T = ±1 and A^ odd/even for some specific 
choices of N. In Figs. and [HI we have given the transition probabilities from Eqs. (|35p - 
(|38p . where for presentation purposes, we have chosen 1 — Pff on the ordinate. From the 
validity requirements in Eq. (1271) we know that Rayleigh-Schrodinger perturbation theory 
will break down at some point. In Figs. [7] and [8], we have therefore presented the curves from 
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numerical calculations. Nevertheless, at least qualitatively, the neutrino transition survival 
probabilities show similar patterns for the analytical and the numerical calculations. 

In what follows, our choice of parameters would correspond in the ADD scenario to a 
4D Planck scale M Pl = 3.4 • 10 18 GeV, a SM Higgs doublet VEV (H) = 174 GeV, a Yukawa 
coupling y — 1 between the active and the bulk neutrinos, and a compactification radius of 
FT 1 = 0.1 eV. In Figs. [5] and El the associated fundamental scale would be M* = 1 TeV, 
which gives mjj ~ 5- lCT 5 eV. In Figs. [7] and [H the corresponding fundamental scale would be 
M* = 50 TeV, which gives mo — 2.5 • 10~ 3 eV. We have distinguished the cases m v = and 
m u ^ 0. In Figs. Eland [HI we have set m v = 0.01 eV, whereas in the other figures, we have 
set m v = 0. We have also made a comparison for the cases T = 1 with the corresponding 
survival probability in the case of a continuous large extra dimension. 

The qualitative behavior of the curves can be understood by looking at the expressions 
for the survival probabilities, i.e., Eqs. (I35l)- fl38l) . From these expressions we observe that 
the dominant effect will be given by the lowest lying modes. For practical purposes one 
can then average over the higher modes, which means that the essential behavior will be 
determined by only a few low modes. Let us first consider the case m v = 0. We obtain 
for T = — 1 and N even the survival probability, when only the first mode is non-averaged, 
as Pff = 2e 2 Bf/K A cos0i + const., where the amplitude and the frequency are given by 
2t 2 B\/K i and 0i, respectively. For the other cases we obtain similar results. 

First, we observe that the frequencies are proportional to 1/R 2 , i.e., a smaller radius 
gives faster oscillations. We also note that the frequencies differ for the cases T = 1 and 
T = — 1. This is due to the different mass eigenvalues that appear in the phases. Thus, we 
have for the case T = — 1, when only the first mode is non-averaged, that the frequency is 
proportional to iV 2 sin 2 n/2N. However, for the case T = 1 the frequency is proportional 
to iV 2 sin 2 ii /N, which is roughly four times larger than for the case T = — 1. This can be 
directly seen in Fig. [51 The figures also show a dependency of the frequency on N. This is 
because the frequency in for example the case T = — 1 is proportional to N 2 sin 2 n/(2N). 
This function grows rapidly for small N and converges quickly to a fixed value, 7r 2 /4. A 
similar relation holds for the case T — 1. This effect is best visible when N < 0(10). 

Second, the amplitude is proportional to m 2 D R 2 so that a change of these parameters 
significantly affects the amplitude. This can be seen by comparing for example Figs. [5] and 
[71 There is also a difference in amplitude between the cases T = 1 and T = — 1. For the 
case T = — 1 the amplitude is proportional to B\, where B n for T = — 1 is given in Eq. f[30|) . 
However, for the case T = 1 the amplitude is proportional to B2 where B n for T = 1 is given 
in Eq. f[39|) . Since Bi(T — — 1) > B 2 (T = 1) the amplitude will be larger for T = —1. 

Note that the case T = 1 and N = 5 differs in a significant way from the other cases. 
The aperiodic behavior of this curve is due to the large interference effect between the two 
lowest modes, which together give the dominating behavior. This effect can be seen in the 
relation between the factors B2 and B%. As we increase N, the effect of B% will be suppressed 
in comparison with B 2 , so that we for large N obtain a sinusoidal-like behavior. For the case 
T = 1 and N = 4 there is no large interference effect between low-lying modes, since the 
sum in this case only includes one term. For the case T = — 1 the corresponding ratio which 
gives the interference effect is B 2 jB\ . This ratio is smaller than the ratio B% / B 2 for the case 
T — 1. Thus, for T = — 1 we do not observe any significant distortion of the periodicity. 



22 



Let us now consider the effect of a non-zero m v . For T = — 1 we note that there will be an 
effect provided that m v R is large enough. If m 2 > l/R 2 , the frequency will be determined 
by 77V However, for the case we have considered, we have chosen m v = 0.01 eV and 
R^ 1 = 0.1 eV, which means that the frequency will mainly be determined by l/R 2 . Thus, 
for the case T = — 1 there will be no drastic changes, which can be seen when comparing 
the upper rows of Figs. [5] and [61 For T = 1, on the other hand, there will be a significant 
effect from m u . This is obvious from Eqs. (13?)) and (|38|) . where the survival probability 
expressions contain a term ae^ 1 , in which a is approximately given by N/(4n 2 R 2 m 2 ,). If 
m u R is sufficiently small, then this term will be the dominating term. Thus, the survival 
probability will be proportional to cos(m 2 L / "IE) . This is the case in the lower row of Fig. 
Note that this effect is due to the presence of a zero mode, which is absent for T = — 1. 

If N is increased, then the curves obtains a more jagged shape because of the interference 
of a large number of KK modes with different frequencies [18]. However, essential properties 
such as the amplitude and the frequency quickly stabilizes. Finally, we observe in Figs. [5] 
and [7| that the case T = 1 reproduces the continuum case [19,21] as expected. 

We have seen that one could, at least in principle, probe T = ±1 as well as the number 
of lattice sites N through neutrino oscillation experiments. For low N the best probe of iV 
is through the frequency. 

7 Summary and Conclusions 

In this paper, we have considered a model for neutrino oscillations in a deconstructed U(l) 
gauge theory defined on the boundary of a two-dimensional disk. If the masses of the link 
fields connecting the center with the boundary are of the order ~ 1 TeV (10 2 TeV), while 
the link fields on the boundary have masses of the order ~ 10 15 GeV (10 12 GeV), then the 
model generates sub-mm lattice spacings between the sites on the boundary. This allows to 
obtain dynamically a non-gravitational large extra dimension with only a few sites. Here, 
we have motivated the general significance of large (sub-mm) lattice spacings by reviewing 
the strong coupling behavior of gravity in local theory spaces for the example of two discrete 
gravitational extra dimensions. 

We have analyzed the mass and mixing properties of a latticized right-handed neutrino, 
which propagates on the boundary of the disk and may be twisted or untwisted. A discrete 
cyclic symmetry, which can be gauged, allows to treat the latticized right-handed neutrino 
as a Wilson fermion with vanishing bare mass. At the same time, the cyclic symmetry also 
introduces in the non-local theory space a local Yukawa interaction between the active SM 
neutrinos and the latticized right-handed neutrino. As a consequence, the model simulates 
key features of the 5D continuum theory for neutrinos in the ADD scenario. 

We have studied the neutrino oscillation effects of the latticized right-handed neutrino 
in terms of the survival probability Pff of a single active flavor / for the cases of a twisted 
(untwisted) lattice fermion and an even (odd) number of sites on the boundary. By taking 
the continuum limit, we could exactly reproduce known oscillation patterns of existing 5D 
continuum theory models. The most direct probe of our model parameters is through the 
frequency of Pff. For example, if the number of lattice sites N is small, then Pff can exhibit 
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a strongly aperiodic behavior for odd N. Possible "odd-even artifacts", however, quickly 
disappear when N becomes large. Generally, twisted and untwisted field configurations can 
be distinguished through the different associated frequencies of Pff, which is for an untwisted 
neutrino roughly four times larger than for a twisted one. The presence or absence of an 
active Majorana neutrino mass also affects the oscillation patterns of twisted and untwisted 
neutrinos in qualitatively different ways. Generally, it should be noted, however, that in 
more elaborate models one would have to include three flavors (as well as matter effects). It 
could also be necessary to take into account additional large extra dimensions. Therefore, 
the results obtained from our model should be viewed as comparatively qualitative. 

The neutrino oscillation effects that are introduced by the KK neutrinos could, in princi- 
ple, be observed in present and future precision neutrino oscillation experiments, such as for 
example KamLAND [49], Borexino [54], or the proposed Double-CHOOZ [55] experiment. 
Borexino would be capable to search for new solar neutrino oscillation effects in an energy 
range E < 1 MeV not covered by Super-Kamiokande or SNO [47]. Our model could be 
tested at short baselines by (future) V e (or u e ) disappearance experiments with sensitivities 
for mixing angles < 0.2 between the active and the KK-neutrinos. Here, it could prove 
useful to employ also two- reactor-two-detector-setups [56], where one may perform mea- 
surements practically free from the typical systematic uncertainties in the reactor neutrino 
fluxes. More generally, one can consider any experiment, which probes the effect of sterile 
neutrinos, provided that one can identify the masses and mixings properly. 

The non-zero mixing between the SM Higgs H and the scalar link and site variables will 
lead to invisible decays H — > WW (if these processes are kinematically allowed), which can 
be checked at the LHC or a future linear collider. 

It is clear, that standard Big Bang nucleosynthesis [57] will be affected by the presence 
of the KK neutrinos. However, the bounds from measurements of the 4 He abundance can be 
alleviated by assuming a primordial lepton asymmetry [58] or with low reheating temperature 
[59]. The constraints on the effective number of neutrino species from large scale structure 
data in conjunction with cosmic microwave background measurements [60] may also be 
evaded by such a lepton asymmetry [61]. Also note that, in this paper, we have assumed the 
same constraints that apply to continuous large gravitational extra dimensions, but it has 
been argued [2] that several of the standard constraints could be relaxed for non-gravitational 
deconstructed dimensions. 
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Figure 5: The neutrino transition survival probability Pff as a function of L/E for some choices 
of TV even or odd and T = ±1. Here we have set R^ 1 = 0.1 eV, mo = 5 • 10 -5 eV, and m u = eV. 
Upper left panel: T = — 1 and TV odd for TV = 5 (dashed curve) and iV = 55 (dotted curve). Upper 
right panel: T = — 1 and TV even for iV = 4 (dashed curve) and TV = 44 (dotted curve). Lower left 
panel: T = 1 and TV odd for iV = 5 (dashed curve), TV = 55 (dotted curve), and the continuum 
model (solid curve). Lower right panel: T = 1 and TV even for TV = 4 (dashed curve), TV = 44 
(dotted curve), and the continuum model (solid curve). 
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Figure 6: The neutrino transition survival probability Pff as a function of L/E for some choices of 
N even or odd and T = ±1. Here we have set i? -1 = 0.1 eV, tod = 5 • 10~ 5 eV, and m v = 0.01 eV. 
Upper left panel: T = — 1 and N odd for N = 5 (dashed curve) and TV = 55 (dotted curve). Upper 
right panel: T = — 1 and ./V even for N = 4 (dashed curve) and TV = 44 (dotted curve). Lower left 
panel: T = 1 and N odd for iV = 5 (dashed curve) and N = 55 (dotted curve). Lower right panel: 
T = 1 and iV even for N = 4 (dashed curve) and iV = 44 (dotted curve). 
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Figure 7: The neutrino transition survival probability Pjf as a function of L/E for some choices 
of iV even or odd and T = ±1. Here we have set R^ 1 = O.leV, vtld = 2.5- 10 -3 eV, and m u = OeV. 
Upper left panel: T = — 1 and N odd for N = 5 (dashed curve) and TV = 55 (dotted curve). Upper 
right panel: T = — 1 and N even for iV = 4 (dashed curve) and iV = 44 (dotted curve). Lower left 
panel: T = 1 and N odd for iV = 5 (dashed curve), N = 55 (dotted curve), and the continuum 
model (solid curve). Lower right panel: T = 1 and N even for N = 4 (dashed curve), N = 44 
(dotted curve), and the continuum model (solid curve). 
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Figure 8: The neutrino transition survival probability Pff as a function of L/E for some choices of 
N even or odd and T = ±1. Here we have set i? -1 = 0.1 eV, mo = 2.5-lCT 3 eV, and m v = O.OleV. 
Upper left panel: T = — 1 and N odd for N = 5 (dashed curve) and TV = 55 (dotted curve). Upper 
right panel: T = — 1 and N even for N = 4 (dashed curve) and N = 44 (dotted curve). Lower left 
panel: T = 1 and iV odd for N = 5 (dashed curve) and N = 55 (dotted curve) . Lower right panel: 
T = 1 and A/" even for iV = 4 (dashed curve) and N = 44 (dotted curve). 
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Figure 9: Amplification of the lower left panel in Fig. (7J comparing the continuum model (solid 
curve) with iV=55 and T = 1 (dotted curve) 
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Figure 10: Triangle diagrams leading to cubic (i) and gauge-gravitational (ii) anomalies in the 
models (a) and (b). These anomalies cancel between neighboring sites. 

A Cancellation of anomalies 

The discrete Zqm symmetry in Eq. (j!2p introduces on each site on the boundary of the disk 
a parity-asymmetry between v nR and u^ R . The model in Sec. H] for the Wilson fermion is 
therefore chiral. If we wish to gauge the Zqm symmetry, we will have to ensure that the 
model remains free from chiral anomalies and that all anomalous contributions from triangle 
diagrams to the three-gauge-boson vertex functions cancel. In the low-energy effective theory, 
the apparent unbroken U(l) gauge symmetry0of the deconstructed model would be c7 (l)dia g 
and it is only at short distance scales that we become aware of the underlying enlarged 
U(1) N+1 gauge symmetry. It is interesting to compare in these two limiting cases the formal 
cancellation of anomalous diagrams by defining the gauge symmetry (containing the discrete 
factor) of our model to be either 

model (a) : G S m x C/(l) diag x Z 6M or model (b) : G SM x U{1) N+1 x Z 6M . (41) 

Implicitly, model (a) becomes equivalent with a non-linear sigma model approximation of 
the deconstructed model in Sec. El Of course, if we are interested in the UV completing 
linear sigma model description, only the anomalies calculated in model (b) are of relevance. 

To work out the anomaly cancellation for models (a) and (b) in Eq. (jHJ), we will denote 
the U(l) n , t/(l)diag, and Z 6M charges of a field / by <&(/), q d m g {f), and 

Qz 6M {f)j respectively 

First, we observe that the fermions / located on the center of the disk - i.e., the SM fermions 
and the fields N a - satisfy qo(f) = 9diag(.f) = 1z 6M if)- Thus, the SM fermions and the N a 
do not contribute to any anomalies and we can from now on concentrate on the anomalous 
diagrams involving only the right-handed neutrinos v Rn and v c Rn . Since the £7(l)di ag and 
U(1) N+1 gauge bosons couple equally to v Rn and v Rn on each site on the boundary, all 
anomalies which do not involve a Zqm coupling vanish automatically. Consider now the 
triangle diagrams in Fig. [TD], which do not have any U(l)i or {7(l)di ag gauge bosons at their 
vertices. If we choose in Eq. ffl2l) M = (N + 2) 2 — 4, then in both models (a) and (b), the 
cubic [^m] 3 anomaly (i) and the gauge-gravitational anomaly (ii) are proportional to 

N N 

Yl^M^Rnf + Qz 6M ( u Rnf] = and ^[gz 6M (^n) + <?z 6jU KJ] = 0, (42) 

ra=l n=l 

4 We neglect here the symmetry breaking introduced by the scalar site variables on the center (see Sec. [3]). 
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the C/(l)diag[-^6Af] 2 anomaly (ii) vanish when summing over all sites. 



where we have used that qz 6M ( u R(n~i)) = -Qz 6M {^ Rn ) and Qz 6M {^m) = qz 6M ( u RN), implying 
that the anomalies cancel between neighboring sites and vanish when summing over all 2N 
right-handed neutrino species. 

The anomaly cancellations differ substantially between model (a) and (b) when evaluating 
the triangle diagrams of the type shown in Fig. [TT], which have at least one U(l) gauge boson 
at one of their vertices. Let us first restrict to model (a). In Fig. [[TJ the mixed Z eM [U (l)di ag ] 2 
anomaly (i) and the f/(l)dia g [-^6J\/] 2 anomaly (ii) are proportional to 

N N 

Y^fosM (vRn) + qz 6M (i4tn)] = and ^[<?z 6M (v Rn ) 2 - qz 6M KJ 2 ] = 0, ( 43 ) 

n=l n=l 

where we have used that gdiag(^n) = — 9diag(^n) — — 1- Again, the diagrams cancel between 
neighboring sites and vanish when summing over all sites. In total, we thus find that all 
divergent triangle diagrams in model (a) formally add up to zero. Let us next consider the 
corresponding anomalies in model (b). The dangerous mixed Zqm\U{1)ti] 2 and U{\) n [ZQM\ 2 
anomalies are obtained from the diagrams in Fig. [11] by replacing in (i) and (ii) the gauge 
bosons according to t/(l)dia g — ► U{l) n . In this case, the diagrams (i) and (ii) in Fig. [11] 
become divergent and the summation over all sites does not remove the divergences, since 
the diagrams belonging to different lattice sites have different external legs and are thus 
inequivalent. 

In order to remove the anomalies in model (b), we add on each site of the boundary 
extra fermions with appropriate quantum numbers. We place on the site corresponding to 
the gauge group U{l) n {n = 1, 2, . . . , N) three additional Dirac spinors, which are written 
component-wise in the Weyl basis as \I/ n = (u n , v n ) T , X n = (r] n ,Vn) T ^ an d X n = (fj n ,Tj n ) T . 
The fields X n , and X n carry the U(l) n charges +1, — (n+ 1), and n + 2, respectively, and 
are singlets under Gsm and the other gauge groups U{l)i^ n . In addition, we assume that 
the extra fermions carry specific Zqm charges. The U(l) n and Zqm charge assignment for all 
SM singlet fermions is summarized in Table [2] 

^From Table [2] we find that the cubic and gauge-gravitational U{l) n and Z 6M anomalies 
still vanish, since the U(l) n and Zqm symmetries satisfy for the extra fermions relations 
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field 


VRn 






v c 

n 


Vn 


Vn 


fjn 


Vn 


q n 


-1 


+1 


+ 1 


-1 


-(n+1) 


n+1 


n + 2 


-(n + 2) 


Qz 6M 


(n + 2) 2 


-(n + 1) 2 


(n + 2) 2 


-(n + 1) 2 


+1 


+1 


-1 


-1 



Table 2: U(l) n (n = 1,2, .. . ,N) and Zqm charges of the fermions on the boundary of the disk. 
These fields transform trivially under the other gauge groups U(l)i^ n . The Zqm charges have been 
normalized with respect to the SM leptons in Eq. (|12|) . 



similar to Eq. (1421) . The mixed ZQM[U(l) n } 2 anomaly is proportional to 

N 

Yl {[^K) - Vn(Vn) 2 ] + lqZ tiM (v Rn ) + QniVnY] 
n=l 

^VZ, M (Vn)-qn(V n ) 2 ] + \qZ 6 AK)+(ln(Vn?}} = 0, (44) 

where each bracketed term inside the sum is zero. Therefore, these anomalies cancel on each 
site. The mixed U(l) n [Z eM } 2 anomaly is proportional to 

N 

E{[-^ 6 A>to) 2 + qZ 6M (i>n) 2 ] + [q n (Vn) + QniVn)} 
n=l 

+ [qZ 6M (v C Rn) 2 -qZ 6M (Kf} + [qn(Vn)+qn(Vn)}} = 0, (45) 

where each bracketed expression inside the sum vanishes. Again, all anomalies cancel in- 
dividually on each lattice site. In total, the Zqm symmetry is therefore anomaly-free. In 
addition, this model is chiral, since the U(1) N+1 x Z 3 m symmetry forbids any bare mass 
terms for the fermions. 



B Neutrino mixing matrices 



In the basis (z/, Uri, Vr 2 , . . . , vrn, v c rx , v c m , . . . , v c RN ), the total neutrino mass matrix M de- 
scribed by the action density £^ in Eq. ([15]) takes the form 



M 



M L M D 
Ml 



(46) 



where Ml and Mp denote (N + 1) x (N + 1) and (N + 1) x N matrices respectively, which 
are explicitly given by 



Mr 



( m u \/Nm D 
Nm D 




• 



o\ 






Mn = U 



( 
-1 +1 
-1 +1 



T 



-1 +1 
-1/ 



(47) 
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Here, the Majorana-like matrix Ml is defined in the basis (u, vri, Vrh ■ ■ ■ > v rn) 
the Dirac-like matrix M D is spanned by (u, v R1 , u R2 , ■ ■ ■ , vrn) and 



whereas 
rn)- In 



Eq. (1461) . "0" denotes an JV x iV matrix with zero entries only. The neutrino masses and 
mixing angles can be determined from the product MM\ which in this basis explicitly reads 



/ 



MM t 



u 



m%+Nm? D 


\/7V m^mjj 




-1 

2 




-1 






—T 


\fNrri£, 
















T 

u 

\TNmuVnu 

T 




5 

N'tri z 
2 + 1™D 
1 u 1 

-1 




u 





u 












-1 


2 


-1 

















-T 






-1 


2 














vJVmo 

u 

\/Nmn 
u 
















2 


-1 




—T 
















-1 


2 


-1 

























-1 


2 -1 



















— T 






-1 2 



\ 



\ ••• —T -1 2 J 

(48) 

where T = ±1 and the blank entries are zero. Next, we want to diagonalize this matrix. 
Since m U) ^/Nttid <C u, we can define the quantity e = ^/Nttid/u <C 1, which we will use 
as an expansion parameter in perturbation theory. We will diagonalize the matrix MM^ in 
steps. First, we perform the transformation MM' — > I' 1 M M I ' using the block- diagonal 
(2N + 1) x (2iV + 1) matrix U = diag(l, U\, U\), where U\ denotes a unitary N x N matrix. 
For T = — 1 and N even, the matrix U\ reads 



/ 





sin 
sin 



N 
A 



ysin 



(JV-l)ir 
JV 




sin 
sin 



A r 
A 



sin 



3(JV-1)tt 
A? 



sin 
sin 





(JV-l)x 
N 

2{N-1)tt 
N 



sm 



(N-l)(N-l)n 
N 



whereas for T = — 1 and iV odd, we have 



cos 
cos 



1 

(JV+1)tt 
N 

2(N+1)it 
N 



COS 



(N-l)(N+l)ir 
N 



COS 
COS 



1 \ 

2(2W-1)tt 



COS 



(JV-1)(2JV-1)t. 



N 



(49) 



/ 







2tt 
A f 



(JV-I)tt 

AT 





sin 22 
sin^ 



3(JV-1)tt 

AT 



sm 
sin 





(JV-2)x 
JV 

2(JV-2)7r 
N 



(N-1)(N-2)tt 
N 



1 

1 
V2 
1 
V2 



1 

V2 



cos 

COS 



1 

(JV+2)tt 
Af 

2(N+2)-n 
N 



COS 



(JV-1)(JV+2)tt 
JV 



COS 
COS 



1 

(2JV-1)tt 
N 

2(2JV-1)7T 

TV 



(JV-1)(2JV-1) 

(50) 
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Similarly, for T 
(3* 



■J2 
1 

1 



\V2 



1 and A" even, we have 




sinf 



2(jV-l)ir 
JV 



sm 
sin 





(JV-2)tt 
JV 

2(JV-2)tt 
iV 



sm 



(JV-1)(JV-2)7T 

JV 



and finally, for T 



V2 
1 

V2 
1 

V2 



\V5 



1 and AT odd 




sinf 



JV 



sin 



2(JV-l)x 

iV 



we have 



sm 
sin 





(jV-l)TT 

JV 
JV 



1 

1 

1 



1 



COS 



COS 



1 

(JV+2)tt 
JV 

2(JV+2)?r 
JV 



COS 



(JV-1)(JV + 2)7T 

JV 



COS 
COS 



1 

(JV+1)7T 

JV 

2(JV+1)^ 
JV 



sm 



(JV-1)(JV-1)7T 

JV 



cos 



(JV-1)(JV+1)7T 

JV 



cos 



cos 



(2JV~2)7r 
JV 

2(2JV-2)tt 
JV 



\ 



COS 



(JV-l)(2JV-2)7r 



JV 



(51) 



cos 

COS 



1 

(2iV-2)7r 
JV 

2(2JV-2)tt 
JV 



COS 



(JV-1)(2JV-2)7T 
JV 



(52) 

In Sec. we consider the case T = — 1 and iV even. The other cases follow in similar ways. 
The rotation from the matrix U\ in Eq. (|49p has the effect of diagonalizing the "gauge-boson- 
like" submatrices in Eq. (|48[) and leads in the new basis to the matrix in Eq. (fT§l) . which 
can then be further diagonalized using perturbation theory as described in Sec. [5J 
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